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X-ray Hardness Variations as an Internal/External Shock 

Diagnostic 

Nathaniel R. Butler^'^ and Daniel Kocevski^ 
ABSTRACT 

The early, highly time-variable X-ray emission immediately following GRBs 
exhibits strong spectral variations that are unlike the temporally smoother emis- 
sion which dominates after t ~ 10^ s. The ratio of hard channel (1.3-10.0 keV) 
to soft channel (0.3-1.3 keV) counts in the Swift X-ray telescope provides a new 
measure delineating the end time of this emission. We define as the time at 
which this transition takes place and measure for 59 events a range of transition 
times that span 10^ s to 10^ s, on average 5 times longer than the prompt Tgg 
duration observed in the Gamma-ray band. It is very likely that the mysteri- 
ous light curve plateau phase and the later powerlaw temporal evolution, both 
of which typically occur at times greater than Th and hence exhibit very little 
hardness ratio evolution, are both produced by external shocking of the surround- 
ing medium and not by the internal shocks thought responsible for the earlier 
emission. We use the apparent lack of spectral evolution to discriminate against 
proposed models for the plateau phase emission. We favor energy injection sce- 
narios with a roughly linearly increasing input energy versus time for six well 
sampled events with nearly fiat light curves at t ^ 10^ — 10^ s. Also, using the 
transition time Th as the delineation between the GRB and afterglow emission, 
we calculate that the kinetic energy in the afterglow shock is typically a factor 
of 10 lower than that released in the GRB. Three very bright events suggest that 
this presents a missing X-ray flux problem rather than an efficiency problem for 
the conversion of kinetic energy into the GRB. Lack of hardness variations in 
these three events may be due to a very highly relatavistic outflow or due to a 
very dense circumburst medium. There are a handful of rare cases of very late 
time t > 10^ s hardness evolution, which may point to residual central engine 
activity at very late time. 



^Townes Fellow, Space Sciences Laboratory, University of California, Berkeley, CA, 94720-7450, USA 
^Astronomy Department, University of California, 445 Campbell Hall, Berkeley, CA 94720-3411, USA 



- 2 - 



Subject headings: gamma rays: bursts — supernovae: general — X-rays: general 



Introduction 



The X-ray telescope (jBurrows et al.ll2005bl ) on Swift (iGehrels et al.ll2004l ) is opening a 
new window into the early lives of 7-ray Bursts (GRBs) and their afterglows. Although hints 
and probable examples of highly time variable X- ray behavior at early time were seen prior 



to Swift (e.g.. |Piro et al.ll2005l : IWatson et al.ll2006l ). the XRT has shown us that this behavior 



is the norm. Nearly all afterglows show a period of rapid flux decline after th e prompt or 



flare emission and about half show bright X-ray flares (e.g.. iBurrows et al.ll2007l ). How these 



observations are to be reconciled with t he well-tested inte r nal /external shock GRB and 
afterglow model (IRees fc Meszarod Il994j : ISari fc PiranI 119971 : ISari. Piran. fc Narayaru Il998 



Wijers fc Galama 



199% 



powerlaws at late time 



- which explained very well pre-Swift observations of simple fading 
comprises a set of key open questions. 



An accurate acco unting of the GRB an d afterglow phenomenology is critical for com- 
parison to the models. lO'Brien et al.l (120061 ) have shown that the late GRB as measured by 
the Burst Alert Telescope (BAT) transitions in time smoothly into the early X-ray counts 
detected by the XRT, provided a correction is made for the different energy bands. This 
de monstrates a close c onnec tion between the early X-ray emission and that from the GRB. 
In lButler fc Kocevskil (120071 . BK07), we fit the BAT and XRT spectra to explicitly show that 
the best-fit models at early time are those which fit GRB spectra well. The early X-ray spec- 
tra look like GRB spectra (but have uFi, peak energies -Epeak in the X-r ay band rather than 



the 7-ray band) a nd evolve spectrally in a similar fashion (see, also, iFalcone et al.l 12006 



Godet et all l2007|l_ Combined with studies in the time domain indicating f ine timescale 



variability (e.g., 



Burrows et al. 



2005a 



Falcone et al.ll2006l : iRomano et al.ll2006l : iPagani et al. 



20061 : iKocevski. Butler, fc BloomlboOTl T we are becoming confident that the X-ray emission 



prior to about lO^s is due to the GRB. The fiat or "plateau phase" light curve which is 
typically present after this phase remains, however, largely mysterious. 

Several models have been proposed to explain the plateau phase light curve. Be- 
cause it is difficult to produce so fiat a decay in the external shock picture, the energet- 
ics m ay be driven by a r e-injection from the central engine or late time internal shocks 



Panaitescu 



Nousek et al.l I2OO6I : IZhang et al.l I2OO6I : iGhiseUini et al 



200 



Off-axis external shocks (Eichler &: GranotI 



200jJPanaitescu et al.l I2OO6 
2006h 



the reverse shock 



(iGenet. Daigne. &: MochkovitchI l2007l: Uhm fc Beloborodovl l2007l ). or time- varying micro- 
physical parameters ( Granot fc Kumar I2OO6 ) may also be responsible. 
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Observationally, IWillingale et al.l (120061 ) have shown that the prompt and afterglow 



emission can be separated near the start of the plateau by fitting two mo dels (of the same 



form) to each inferred component. This split falls short of decisive because IWillingale et al. 



( 120061 ) are unable to measure the rise of the afterglow component from under the prompt 
component. The number of degrees of freedom in the model is large (6-8) and compara- 
ble to the typical number of powerlaw segments in broken powerlaw fits, which assume no 
separation into prompt and afterglow components. As we show below, a cleaner separation 
that requires no manual removal of flare-like emission is possible if we consider the spec- 
tral variations at early time. This can be demonstrated through the use of time resolved 
spectroscopy as discussed above, although such efforts are limited to bright events with high 
signal-to-noise. Because spectral fits are not required, variations in the X-ray hardness ratio 
provides an alternative. Studying the hardness ratio, we can link the plateau phase emission 
to late-time external shock emission for even faint bursts. 

After a brief review of the X-ray phenomenology versus time gleaned from powerlaw fits 
(Section [3]), we discuss in Section H] how the hardness evolution implies a separation between 
prompt and afterglow emission. Stable hardness ratios during and after the plateau phase 
are exploited to constrain the GRB and afterglow models in Sections O and El 



2. Data Reduction 

We download the Swift XRT data from the Swift Archiv^. The data are processed with 
version 0.10.3 of the xrtpipeline reduction script from the HEAsoft 6.0.(1^ software release. 
We employ the latest (2006-12-19) calibration files. The reduction of XRT data from cleaned 
event lists output by xrtpipeline to science ready light curves and spectra is described in 
detail in Butler & Kocevski (2007). Our final light curves have a fixed signal-to-noise of 3 
in the 0.3-10.0 keV band. 

We define an X-ray hardness ratio HR as the fraction of counts in the 1.3-10.0 keV 
band to the counts in the 0.3-1.3 keV band. On average, this ratio is equal to unit y for 



XRT data. The mean energy index (flux proportional to E~^) is /? = 1 (lButlerll2007al ). We 
show in BK07 that the column densities Nh, as inferred from soft X-ray absorption, do not 
appear to change in time for Swift afterglows. To lowest order for a typical column density 
Nh = 10^1 cm-2, (3^1- 0.9\og^{HR). 



^ftp: / /legacy.gsfc. nasa.gov/swift /data 

^http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
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2.1. Light Curve Region Selection and Fitting 

To group the data into separate regions of similar temporal and spectral evolution, we 



fit the data using an extension of the Bayesian blocks algorithm (IScargldll998l ) to piecewise 
logarithmic data. Our implementation is simple and requires no human intervention. Con- 
sidering each data point as the location of a possible powerlaw break in the light curve, we 
calculate xl fo^' every possible connection of the data points. This search can be done effi- 
ciently using publicly available Bayesian blocks cod^. Because we also include an additional 
term (a prior against the break) with each possible new segment of Axl = 9.2, each new 
segment must improve the fit at ^ 99% confidence. The Axl additions exclude models with 
many breaks, while the final fit without the Axl contributions has xl ^ ^ typically 3-5 
powerlaw segments. 

So that noise fiuctuations in the data do not generate spurious lig ht curve breaks, we 



denoise the soft and hard channel light curve with Haar wavelets (see, iKolaczyk fc Dixon 



20001 ) prior to fitting. We fit the count rate and hardness simultaneously so that spectrally 
distinct regions are separated. The fits are plotted in red in the figures below. Light curves 
and fits for all Swift events can be downloaded fron]^. 



3. Temporal/Spectral Overview of Swift and pre-Swift X-ray Afterglows 

The X-ray spectral and temporal properties of GRB afterglows are typically garnered 
from powerlaw fits in time and energy. For Swift data, this can be done for individual events 
at multiple epochs. A detailed study of the evolution across all of these snapshots is beyond 
the scope of this paper, and we restrict our analysis in this section to the time regions as 
though they were separate events. It is reasonable to compare Swift data to pre-Swift data 
in this fashion, because pre-Swift fits correspond to a narrow time windows in the life of 
each afterglow and do not typically allow for time evolution studies for individual bursts. 
We determine fit regions with the blocking algorithm described above, and the results are 
binned in time into the three time epochs. 

We show in BK07 that powerlaw energy fits are inappropriate at early time. This 
can also be seen in the first panel of Figure [1], which plots the energy index /3 versus the 
temporal index a for 59 Swift afterglows prior to and including GRB 070208 (see Table 1). 
For t < lO^s, there is a very wide range in /3, consistent with the combined range of low- 



3 http :/ /space . mit . edu/ CXC/ analysis/ SITAR/ 
"'http: / / astro. berkeley.edu/^nat / swift 
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energ y and high-energy indices observed in BATSE GRBs (jPreece et aLll200Cll : iKaneko et al. 
20061 ). Nearly half of these points (42%) are inconsistent with any of the external shock 
synchrotron models plotted as diagonal lines. The fraction is greater than a half (52%) if 
we consider bursts separately, rather than plotting multiple spect ra from indi v idual events. 
We show in BK07 that the X-ray spectra are well fit by the same iBand et al.l (119931 ) model 
which fits the GRB spectra. The large scatter in the time indices for t < lO^s is due to rapid 
light curve decays and flaring. In Section H] we show that the X-ray hardness can be used to 
infer the end of this phase. 

From 10^ ^t<^ lO^s, the X-ray light curves typically decay at a much slower rate. As 
shown in Figured] (middle panel), there is apparently little spectral evolution. Most of the 
fits (79%) and about half of the total number of bursts (53%) are c onsist ent with an adiabatic 
shock observed above the cooling frequency Uc- IWillingale et al.l (120061 ) also find that about 
half of all bursts are not consistent with this synchrotron model due to anomalously slow 
time decays. The fraction of consistent spectra is larger, because the events without plateaus 
are generally brighter at this stage. The plateau events can be modelled assum i ng a smooth 



re-ini ection of energy into the external shock at late time (iNousek et al.ll2006l : IZhang et al. 



20061) or by decele ration of the external shock in a wind density {n oc R ^) external medium 
dPanaitesculboOTh . 



After lO^s, the fits exhibit a tight clustering in both a and (3 and 88% of the tempo- 
ral/spectral snapshots (or 80% of bursts) are consistent with the behavior expected from an 
adiabatically expanding shock in the circumburst medium, emitting synchrotron radiatio n 



above the cooling frequency with electron index p pa 2 fe.g..lSari. Piran. fc Narayarull998l ). 



The pie- Swift X-ray d ata from Beppo-SAX (see, e .g. Jde Pasquale et al.ll2006l ) and Chandra 
and XMM (see, e.g., iGendre. Corsi. fc Pird l2006l ) are all taken beginning after this time 
and show closely con sistent behavior with th e S wift events p l otted here. Apparent here, but 
discussed in detail in IWiUingale eTaP toodl and iPanaitescul (120071 ). few of these events here 
fit the expectation for a jetted afterglow jRhoadslll999[ ). 



4. T ~ lO^s and the End of Hardness Variations 

Figure [2] plots the X-ray hardness ratio versus time since the GRB trigger for the full 
sample of XRT afterglows. We also show ±1 times the root-mean-square scatter in the data 
in red. The scatter is several times greater prior to lO^s than after. The reason for this 
scatter becomes clear when looking at the afterglows from individual bursts. Figure [3] shows 
the X-ray light curves and coincident X-ray hardness ratios for six events. Each shows an 
early period of strong hardness variation, which flattens out to a late-time value after several 



- 6 - 



hundred seconds. The hard, late-time component appears to overwhelm the soft, early-time 
component. We mark as the Th the time at which the hardness reaches a minimum in each 
plot, before gradually increasing to a constant late-time value. 

Figure H] plots the distribution of end times of these hardness variations Th for 50 
afterglows (also Table 1). The flux prior to Th is typically far softer than that after the start 
of the plateau phase, and the plateau phase becomes evident in HR plots prior to becoming 
evident in the 0.3-10.0 keV count rate. We note that the 0.3-10.0 keV light curves transition 
smoothly across Th- Therefore the energy integrated light curve cannot be used to measure 
Th- Several additional events, with and without flaring light curves, are plotted in BK07. 

There is a significant anti-correlation between the time since the end of the BAT emission 
Th — Tgo and the X-ray flux at Th (Kendall's tk = —0.46, signif.= 2 x 10~^), which reflects 
the decay of the GRB flux and its intersection at Th with a range of possible flux levels for 
the afterglow plateau phase component. 

The end times Th are on average 5 ±3 longer than the prompt Tgg durations. The X-ray 
flux at Th is on average io^-^='=''-^ times fainter than the average flux for the GRB measured in 
BAT. These quantities are only weakly correlated: Tgo versus Th has a Kendall's tk = 0.20, 
signif.= 0.04 and Fx versus has a Kendall's tk = 0.03, signif.= 0.71. 

We can use the light cu rve taxonomy devel oped by other authors to re late our times 
Tf{ to the "canonical" (e.g ., iNousek et al.ll2006l ) light curve decay phases. lO'Brien et al. 



hood ) and lWillingale et all J2006h divide the light curves into prompt and afterglow phases 



by fitting models to the energy integrated light curves. These authors estimate durations 
which represent the brightest and most slowly decaying time regions, excluding the rapidly 
time-decaying tails of the em ission episodes. Our Th values are on average ten times longer 
than the prompt time Tp in IWillingale et al.l (120061 ). but there is no significant correlation 
[tk = —0.04, signif .= 0.75, for N = 35 b ursts). There is also no significant correlation 
between Th and the Iwilhngale et all J2006h Ta {tr = 0.1, signif.= 0.43, for = 35), which 
approximately measures the end of the plateau phase and is on average twenty times larger 
than our Th- 

Given Th as a dividing line between emission with strikingly different temporal and 
spectral characteristics — which we can interpret as a dividing line in time between GRB 
and afterglow emission — it is possible to separate and compare the fluence from the GRB 
and afterglow. Figure plots the fluence prior to Th versus the fluence after Th- The 
afterglow emits an amount of energy proportional to {tr = 0.35, signif.= 3 x 10~^) and 
]^qi.o±o.5 ^jjjies lower than the GRB. The fluence in the X-ray band prior to Th contributes 
only 10% additional fluence, on average, to the GRB as observed in BAT. These quantities 
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are possibly weakly correlated {tk = 0.18, signif.= 0.07). The prompt and afterglow fiuences 
we find here are consistent with those reported by Willingale et al. (j2006 . Fig 3) using two- 
component model fits to the light curves, but our fiuences correlate with less scatter (Figure 

ED- 



5. Achromaticity of the X-ray Plateau Phase at 10'^ ^ lO^s 

The events selected for plotting in Figure [3] have prominent flat X-ray light curves at 
10^ <s t <^ lO^s. On the HR panel for each burst in Figure [3] we print the maximal HR 
variation between data points after lO^s. This is always less than 0.5, corresponding to 
SP < 0.5. 60 = 0.5 is the expectation for changes d ue to cooling in the external shock (e.g.. 



Sari. Piran. fc NarayanI Il998l : I Chevalier fc Lil |2000| ) . If we place a tighter constraint on the 
allowed hardness ratio variations by allowing only powerlaw increases or decreases after lO^s, 
then the limits on 6f3 are much tighter, 6(3 ^ 0.1. We thus see no evidence for significant 
spectral evolution in these events after lO'^s. 

We also observe no significant variation in j3 from powerlaw fits, consistent with the HR 
analysis. Combining this information to derive a constant (3 in time for each event, we can 
plot a and (3 across the break (Figure E]). Many of the light curve breaks in Figure [3] are 
gradual and are fit here with multiple powerlaw segments. As seen in Figure Eland discussed 
more below, the values of (a,/3) during the plateau are consistent with values produced in 
energy injection models. As the light curve breaks, the {a, (3) approach those expected from 
external shock models without energy injection. 



6. Discussion 



We have exploited an autonomous spectral/temporal division of early afterglow data 
to isolate the time when models for the GRB afterglows well-tested prior to the launch of 
Swift first begin to match the data well. From powerlaw fits in time and energy (Section [3]), 
the afterglow model s appear to break down strongly prior to t ^ W^s. This is consistent 
with the findings of lO'Brien et al.l (120061 ): IWillingale et al.l (120061 ). We focus here on rapid 
time variations in the X-ray hardness, which end by Th ~ 10^ — 10^ s and therefore allow 
for a clean separation of early, GRB-like emission and later afterglow-like emission without 
hardness variations. 



Our finding here that the end time of hardness variations anti-correlates strongly 
with the X-ray flux at Th (Figure H]) likely has a trivial explanation: because the light curve 
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is decreasing logarithmically after Tgo, the length of the duration Th — Tgo simply reflects 
the faintness of the afterglow. More interesting, the large dynamic range in this correlation 
between prompt and afterglow fiuences (Figure E]) indicates that some physical feature of the 
explosion or circumburst cite must be able to substantially modulate the fraction of energy 
in highly relativistic material (the GRB) or the shock kinetic energy (the X-ray afterglow). 

GRB models must be able to explain how the fireball deceleration (see, Meszaro^ 2002 



and references therein) can be postponed until after t ^ 10^ s — and probably until after 
t ^ 10^ s — when the afterglow light curve is no longer fiat or rising. There must be no 
apparent imprint of the external medium on the light curve or spectrum prio r to t ~ 10 ^ 



s. Because the deceleration time tr\p. r = 390(1 + ^)-£^shock sa'^i ^^^-'^2 ^''^ s (e.g., iPiranI 11999 



Ramirez- Ruiz. Merloni. Sz Reesll200ll ). this might be accomplished by increasing in the shock 



kinetic energy -Eshock.ss [10^^ erg] via energy injection, by having a very low density n [cm~^], 
or by placing much of the outflow in low Lorentz factor F = 100 F2 material. The importance 
of these parameters becomes more apparent if we focus on extreme cases (Sections 16.11 and 
16. 2p or on modelling of the more typical cases (Section l6.3p . 



6.1. Events with No Hardness Variations 



There was an expect ation, based primar i ly on detections of putative afterglow compo- 

nents in the tails of GRBs ( Connaughton 2002; Giblin et al.ll2002 ; Lazzati. Ramirez-Ruiz, fc Ghisellini 
2OOII ) and also from studies extrapolating the X-ray afterglow flu x back to the prompt emis- 



sion (e.g.. lCosta et al.lll997l : iPiro et al.lll998l : iFrontera et al.ll2000l ). that the X-ray afterglows 
measured by Swift would directly proceed and connect to the prompt emission. That is, the 
emission at t ~ 10^ — lO^s was expected to mirror the temporal/spectral properties of the 
well-established later emission plotted in the third panel of Figure [H Instead, the rapid time- 
variation and plateau phases generally occur, with strong spectral variations during the rapid 
time-variation phase. It is a well known fact that the late-time emission in Swift GRB X-ray 
after glows does not typically extrapolate back in time to the end of the prompt emi ssion 



le.g- 



Nousek et al.ll2006l : IZhang et al.ll2006l : lO'Brien et al.ll2006l : IWillineale et al.ll2006r ). In 

stead, the X-ray light curve is relatively flat prior to t ~ lO'^s, and this resul ts in typical 



X-ray f iuences that are ten tim es lower than the GRB fluence (Figure [5], also, IZhang et al. 



2006bl : iGranot fc Kumad 120061 ) . Bec ause the late time X-ray emission is typica lly argued 
to trace the kinetic energy explosion (jKumarll2000l : iFreedman fc Waxmann II2OOII ). this may 
lead to an efficie ncy problem because the shock energy available to the GRB is very low 
(e.g., 



PiranI 119991 ). 



Because efficiency of the conversion of shock kinetic energy to the energy released by 
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the GRB is an intrinic property, we would expect that the efficiency should vary little from 
event to event and that all events should exhibit an X-ray plateau. However, there are few 
bursts — <10% of the 30 or so afterglow observations which began early and measured 
with high signal-to-noise the tail of the prompt emission — which do not show early HR 
variation and do appear to show an afterglow-like component decaying as a p owerlaw after 
the p rompt emissio n at t >, 100s. T he best example GRBs 050717 (also. 



2006 ). 060105 falso. lTashiro et al.ll2006h . and 061007 falso. lMundell et al. 



2007 



Krimm et al. 



Schady et al. 



20071 ) are shown in Figure El These have especially bright and hard prompt emission with 
-Epeak ^ 500 keV (suggestive of high V) and hard X-ray emission detected beginning after 
t 90s with energy index /5 ~ 1. Perhaps an early deceleration occurs for these events due to 
an anomalously high circumburst density. Could these events also be telling us that energy 
is present but not observed in the soft X-ray band in the other, more common afterglows 
with prominent plateau phases? Where is this energy? 



One possibility is the two jet model (lEichler fc Granot 



20061). with a GR B jet containing 



more kinetic energy per solid angle than in the afterglow. iPanaitescul (120071 ) has tested and 
found no evidence for this scenario. Another possibility is that the early forward shock 
emission is suppressed as i t scatters (to the GeV-TeV range) photons from the late-time 
internal shocks and flares (IWang. Li. fc Meszarosll2006l ). Understanding a possible energy 
removal mechanism and its impact on the fluence correlation (Figure ISJ will likely help also 
to understand a p ossible correlation be twee n the end of the plateau phase and the GRB 
energy reported in lO'Brien et al.l (120061 ) and iNava et al.l (120071 ). 



Finally, we note that GRBs 050717, 060105, and 061007 may be representative of a 
separate class of GRBs to which previous missions were more sensitive. This would explain 
why afterglow-like tails are rarely observed just following Swift GRBs. In part, previous 
mission were likely also less sensitive to the very soft emission detected by Swift. A handful 
of these — 050502B, 050724, 061222B, 070129 — have very soft emission between Tgo and 
Th which is gr eater than the pro mpt fluence. These events may help us to understand how 
X-ray Flashes (iHeise et al.ll2000l ) are related to GRBs. A late and bright X- ray flare as in 



050502B likely produced the soft X-ray excess in the enigmatic GRB 031203 (IWatson et al. 
2006h . 



6.2. Events with Very Late-time Hardness Variations, The External Shock? 

There are rare events which show hardness variations with AHR > 0.5 after t ^ lO^s 
(Figure 121), which is probably too late for an explanation involving the deceleration of the 
GRB ffieball. The unusual supernova-GRB 060218 stops varying in hardness just after lO^s. 
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(Give n the extremely long prompt duration of this event (T90 ~ 2 x lO^s; ISakamoto et al. 
20061 ). the late hardness evolution may not be unusual. This GRB and afterglow produces a 
clear arc from t ^ 300 — 3000s in Figure [21 which demonstrates a HR values and evolution 
distinct from those observed in any other event.) The high-z GRB 050904 light curve appears 
to consist entirely of flares, and this is reflected in late-time hardness variations. 

The hardness increases at late-time for GRB 060206 (the outlier in Figure H]). The 
additional cases (GRBs 050315, 060105, 060814) show a decreasing hardness on a timescale 
similar to the observations time. The hardness variation is consistent with the factor of two 
expected fror n the AP = 0.5 change expect e d from a cooling brea k in the synchrotron shock 



picture (e.g., ISari. Piran. fc NarayanI Il998l : IChevalier &: Lill2000l ). Although we believe we 



have accounted for the flux of a nearby source, it is possible that the hardness increase in 
the case of GRB 060206 is due to contamination by that source. This is not an issue for 
GRBs 050315, 060105, and 060814. 



For GRB 060206 at z = 4.045 (IFynbo et al.ll2006[ ). the source-frame GRB energy release 
is i^iso = 4.2lQg X 10^^ erg. The break to increased hardness in Figure [7] for a wind medium 
implies a reasonable density A^, = 0.13(eB/0.01)^^/''. If the density were uniform, as would 
be inferred from the breaks to softer spectra in the other 3 events, a very low density 
n = 1.7x 10~^(eB/0.01)~^''^ is implied. Unless the other 3 events are at low redshift {z ^ 0.1), 
which is unlikely given the lack of bright optical emission in each case, the implied densities 
are anomalously low. 

Falling back on our basic ignorance of the nature of the prompt engine and its timescales 
for energy output, the late time variation could also be due to the central engine. The light 
curves do appear more structured than simple powerlaws. In Figure [71 we fit the observed 
0.3-10 keV rate model to the HR in order to derive a hardness intensity correlation index. 
These are printed in the figure and are typically ~ 0.5 for the early, GRB-like emission and 
^ for the late emission. As we discuss in BK07, a probable explanation for this correlation 
at early times is a relativistic viewing effect due to photon arrival from emitting regions 
off the line of site. The long emission timescales could be due to curvature radiation at 
large radii. It is interesting to note that the late hardness evolution in GR B 050315 rules 



out as sociating the late break with an achromatic jet break, as was done in IVaughan et al. 



6.3. Energy Injection Scenarios 



The lack of spectral evolution during the light curve plateau phase and at later time 
suggests strongly that these episodes are generated by the same emission mechanism. In 
the fireball mo del, this points to a smooth late-time energy injection that refreshes the 
external shock (INousek et al.l l2006l : IZhang et al.ll2006l ). which alters only the time decay 
rate and not the spectral regime in which the X-ray synchrotron emission occurs. The 
external shock without energ y injection cannot produce the observed fiat light curves (e.g., 
Sari. Piran. &: Narayanlll998l ). 



Energy injection i nto t he GRB external shock was first discussed by iPaczyhskil (119981 ) 
an dlRees &: Meszarod (119981) . According to the very general energy injection models outlined 



m 



Nousek et al. 



(120061 ) . the X-ray light curve is propped up by an insertion of energy E in 
time as E o(.t°'. For data observed above the synchrotron cooling break Uc and peak frequency 
Um, the flux varies as Fx oc t~"'"j, with ainj = a — a{(3 + l)/2 = (3 — a)(3/2 — (1 -|- a)/2. 
Curves for a in the range 0.5 — 1.1 — which is needed to account for the well-sampled plateau 
phase events described in Section |5] — are shown in Figure [61 Time flows from the left to 
the right in this figure as the light curve smoothly breaks. We consider the left-most points 
for each burst as those most likely to reflect the true energy injection profile. If the X-ray 
band is below the cooling break, the slopes in Figure E] remain the same for the same a but 
the offsets shift. This leads to acceptable fits with a = 0.7 — 1.1 for a constant density (ISM) 
medium and a = 1.0 — 1.4 for a wind density medium. An E oc t scaling implies a central 
engine with an approximately constant late-time luminosity or an ordered flow of internal 
shock mat erial with M(> T) oc F"^ — F"^ (ISM med ium) reaching the afterglow shock at 



late time (INousek et al.ll2006l : iGranot &: Kumarll2006l ). 



After the break, the fits become consistent with external shock models ranging between 
the expectation for a spherical expansion observed at z/ > i/c, T-'m and a jetted expansion 
in the same regime. This may imply that gradual jet breaks are present in these events, 
although only one (060614) ever reaches the expected late-time decay rate. 

Alternatively, the X-ray band could be in the Vc > v > Vm regime and the energy 
injection could be turning off gradually. The steep decay in the case of 060614 at late time 
requires a wind density medium. However, this can be ruled out from optical data reported 
to the GC^. During the plateau, the optical light curve rise as tO-38±o.o6^ compared to 
^-o.o3±o.o5 measured in the X-ray band. This behavior is consistent with v^. between the 
optical and X-ray band and energy injection with a = 0.69 ± 0.08. There are no optical 



^http://gcn. gsfc.nasa.gov 
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points after the candidate jet break to verify achromaticity. GRB 060729 has consistent 
optical (f-^-'^'^^o-o^^ and X-ray (^-o-26±o.o4^ hght curve indices during the plateau, and both 
light curves break to consistent decays thereafter. The post break decay is consistent with 
expansion into an ISM medium without a jet break. The energy injection prior to the break 
is fit by a = 0.92 ± 0.05 {u > u^, u^). 



Panaitescul (120071 ) propose a very simple model to explain some plateau light curves 



and spectra. For a wind density medium, t^ec is a strong function of the bulk Lorenz factor, 
^dec = 6(1 + z)Es\ioc]i,^2,^~^V2^ s for a typical Wolf-Rayet wind density of 5 x lO^^A^, g cm~^. 
The afterglow will not peak until 10^ — lO^s if, after the internal shocks are through, the 
effective V ~ 20. During dece leration, the flow coasts and the hght curve stays relatively flat, 
a = [3 — 1 (jPanaitescul 120071 ) . for v > Uc, ^m- Contrarily, deceleration by a uniform density 
medium produces a sharply rising light curve, which is not observed. This a,j3 relation for 
the wind medium takes the same form as that for the a = 1 energy injection model, and 
it appears to be roughly consistent with most of the plateaus in Figure [61 In this picture, 
the jet break will coincide with the end of the plateau if the opening angle is ^ ~ l/F ^ 3°, 
which may well occur for some events. 

Several additional models have been proposed to explain the X-ray plateau phase, and 
many of these can be constrained by a constant X-ray spectral slope and from the fact that a 
distinct hardness evolution separates the plateau phase emission from the GRB-like emission 
prior. One possibility — which we can rule from lack of hardness variations because it 
requires contribution from the spectrally varying tail of the GRB — is that off-a xis afterglow 



plus late GRB emission combine to produce the plateau (lEichler fc Granot 



2006h . In a similar 



fashion, we may be able to rule out the "late prompt" model of iGhisellini et al.l (120071 ). 
although spectral variations would be modest in that rnodel due to low F. Models i nvolving 
the reverse shock (IGenet. Daigne. fc MochkovitchI 120071 : lUhm fc Beloborodovl 120071 ) appear 
to produce spectra variations in the X-ray band. Invers e Compton models wh i ch ext ract 
afterglow flux should also change the spectrum, but see IWang. Li. &: Meszarod (120061 ). A 
final, more exotic possibility w hich we cannot rule out i s time evolution of the microphysical 
parameters defining the shock (iGranot &: Kumarll2006l ). 



7. Conclusions 

We have shown that GRB and early afterglow light curve prior to Th ~ lO^s is highly 
time and energy variable. The flux at the end of these variations and during the X-ray 
plateau phase exhibits non-changing X-ray hardness like the latest X-ray afterglow emission 
that is well modelled by a synchrotron external shock [t ^ lO^s; Section [3]). The afterglow 
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flux is typically ten times lower than would be estimated from a simple extrapolation of the 
GRB flux after Tgo ~ Tf//5. Explaining how the GRB deceleration can be postponed until 
after t ^ lO^s is a central challenge to those modelling GRB and their afterglows. 

Observations prior to Swift which imply a common early onset of the afterglow may 
point to a class of bursts rarely observed by Swift. The three Swift examples discussed 
(GRBs 050717, 060105, and 061007) have energetic and hard prompt emission. These bursts 
may be those most rapidly decelerated by the circumburst medium. These events likely 
have the most high energy photons for GLAST to observe, unless photons missing from the 
early afterglows of softer events are preferentiall y up-scattered to high energies by late time 



shocking (see, e.g.. IWang. Li. fc Meszarosll2006l ). In either cases, GLy45T observations will 



be crucial for understanding this diversity and for helping us to u nderstand th e Swift phe- 



nomenology relative to that observed in previous missions (see, e.g.. IZhang||2007l ) . Additional 
long wavelength observations (e.g., in the optical/IR) are also essential at times t ^ lO^s, 
because these better probe the circumburst density structure. 
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Table 1: X-ray Hardness Evolution End Times Th 



GRB Th [s] GRB Th [s] GRB Th [s] GRB Th [s] 



050219A 


165^35 


050315 


991+100 


050502B 


974^1^ 


050607 


410 ±40 


050712 




050713A 


203^i^i 


050714B 


400t^gon 


050716 




050721 


1135_242 


050724 


090+38 


050730 


7791}} 


050801 


707+142 

'•J '-142 


050814 


562_274 


050820A 


4750112^ 


050822 


600+i^2 


050922B 


,-,91 + 129 
iizi_-^29 


051117A 


1659t^^ 


051227 


220 ± 100 


060115 


594t}^^ 


060204B 


378+i^7 


060206 


21284tl;^^^^ 


060210 


426+^ 


060211A 


391+^? 


060312 


1761^^ 


060418 


17511^2 


060427 


2261^6 


060428A 




060428B 


271135 


060510B 


415t|^^ 


060526 


3741^1 


060604 


193+1 


060607A 


2891^ 


060614 


2588_2i27 


060707 




060708 


88t}} 


060714 


232+}2 


060719 




060729 


971 +26 
^' J^-26 


060904A 


413t47 


060904B 


1948t}^^l 


060929 


630;ii 


061110A 


390+^6 


061121 


142+}} 


061202 


390t}^^ 


061222A 


203+1^ 


061222B 


275il« 


070107 


477ii 


070110 


2541}^; 


070129 




070208 


25851^1^^ 






Limits 
















050401 


<153 


050416A 


<126 


051016B 


<91 


051109A 


<195 


060111B 


<145 


060211B 


<95 


060219 


<158 


060502A 


<88 


060510A 


<110 


060906 


<276 


060908 


<90 







Notes: Several afterglows allow for only a Th limit measurement. 
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T<1x10's 10^s<T<2x10'*s T>2x10''s 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 

Time Decay Index a [F « t"] a a 



Fig. 1. — Powerlaw energy index versus powerlaw light curve index for the X-ray afterglows 
at 3-time epochs after the GRB. The solid curves show the expected " closure relations" for an 
adiabatic expanding external shock emitting synchrotron radiation (jSari. Piran. fc Narayan 
19981 : IChevalier fcLjlbnod : ISari. Piran. fc Halpernl[l999l ). Only after t = 2 x lO^s — a time 
corresponding to all pre-Swift X-ray afterglow observations — are the data overwhelmingly 
consistent with the models. 
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Fig. 2. — The scatter in hardness ratio for the full Swift sample prior to GRB 070208 drops 
rapidly with observation time since the GRB trigger. We have subtracted off a constant HR 
from each burst, which we obtain by fitting the data at t > lO^s. 
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Fig. 3. — X-ray Light curve and hardness ratio HR plots for several GRB afterglows 
with prominent light curve "plateau" phases. The temporal indices of the best fit temporal 
powerlaw models (red curves) for the counts rates are shown. Th marks the end of an 
early phase of strong hardness evolution in each case. At later time, HR is consistent with 
constant. The dotted lines mark the best-fit late-time value as well as e°'^ and e~°'^ times the 
best-fit value, as a characteristic expected range for variations in external shock synchrotron 
models. A limit on the maximum deviation (1-sigma) from the best-fit late-time HR value 
is also given. 
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Fig. 4. — Frequency distribution (black histogram) of the latest time at which spectral 
hardness variations are observed Th (Table 1) for 50 events. The hardness typically remains 
constant after this time (Figure [3]). The red histogram for fewer events is divided through 
by the event redshift (1 + z). The subpanel shows the Tgo duration measured in BAT versus 
the hardness variation end-time and also the X-ray flux at Th versus Th- There is a weak 
correlation with Tgo and a strong correlation with Fx- On average, Th is 5 times greater 
than Tgo and the X-ray flux at Th is 175 times below the average GRB flux observed by 
BAT. 
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Fig. 5. — Time integrated flux before Tjy compared to after Th". The fluence in the afterglows 
(T > Th) component is on average ten times less than the fluence in the prompt GRB and 
early X-ray emission. 
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Fig. 6. — Energy and time indices for powerlaw segments in the six bright plateau events 
from Figure [3l Separate events are color coded. Indices for the flat, plateau phase of the 
light curve in each case are consistent with energy injection models (dotted lines shown are 
for V > Uc, I'm)- The break in time to a steeper fade — which occurs roughly after the model 
marked a = 0.5 in this plot — is gradual. The powerlaw indices at late time approach those 
expected from the external shock models. Models for the late time emission of spherical and 
jetted shocks are plotted. 
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Fig. 7. — X-ray Light curve and hardness ratio HR plots for rare GRB afterglows which 
show late-time hardness variations. These variations may imply very late time central engine 
activity. Fits of the observe d flux to HR are plotte d in red, and index values ~ 0.5 are 
expected for internal shocks (iButler fc Kocevskil 120071 ). 
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Fig. 8. — X-ray Light curve and hardness ratio HR plots for three events where the GRB tail 
(blue points) connects directly with a powerlaw (afterglow-like) X-ray light curve exhibiting 
little spectral variation. The X-ray data in these cases are well fit by powerlaws with energy 
index (3 = 1.0 ±0.1 throughout. The dotted lines mark an expected range for variations in 
external shock synchrotron models (also Figure [3]). In the case of GRB 060105, the 1-sigma 
limit on Alog{HR) oc A/3/0.9 is consistent with a possible cooling break at t ^ lO^s. 



